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O . Abstract 

o 

I We simulate the distribution of neutral hydrogen (HI) at the redshifts 2 = 1.3 

and 3.4 using a cosmological N-body simulation along with a prescription for 
2 ' assigning HI masses to the particles. The HI is distributed in clouds whose 

^ ' properties are consistent with those of the damped Lyman-a absorption systems 

(DLAs) seen in quasar spectra. The clustering properties of these clouds are 
, identical to those of the dark matter. We use this to simulate the redshifted 

^ . HI emission expected at 610 MHz and 325 MHz, two of the observing bands a 

the GMRT. These are used to predict the correlations expected between the 
complex visibilities measured at different baselines and frequencies in radio- 
interferometric observations with the GMRT. The visibility correlations directly 
probe the power spectrum of HI fluctuations at the epoch when the HI emission 
originated, and this holds the possibility of using HI observations to study large- 
scale structures at high z. 
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1 Introduction 



Observations of Lyman-a absorption lines seen in quasar spectra are an impor- 
tant probe of the distribution of neutral hydrogen (HI) at high redsliifts. These 
observations show that the bulk of the neutral gas in the redshift range 1 < z < 
3.5 is in HI elouds with eolumn densities greater than 2 x 10^*^ atoms/cm^ (Peroux 
et al. 2001, Storrie-Lombardi, McMahon, Irwin 1996, Lanzetta, Wolfe, & Turn- 
shek 1995). The damped Lyman-a absorption lines produced by these clouds 
indicate U,gas{z)^ the comoving density of neutral gas expressed as a fraction 
of the present critical density, to be nearly constant at a value ^gas{z) ~ 10^^ 
(Peroux et al. 2001). 

In this paper we simulate the HI emission expected from these clouds. The 
aim of the exercise is to investigating the possibility of detecting the redshifted 
HI emission using the Giant Meterwave Radio Telescope (GMRT; Swarup et 
al., 1991). We focus on two of the GMRT frequency bands centered at 610 MHz 
and 325 MHz corresponding to HI emission from redshifts z = 1.3 and 3.4. 
The HI flux from individual clouds (< 10/xJy) is too weak to be detected by 
GMRT unless the image of the cloud is significantly magnified by an intervening 
cluster gravitational lens (Saini, Bharadwaj and Sethi, 2001). Although we may 
not be able to detect individual clouds, the redshifted HI emission from the 
distribution of clouds will appear as a background radiation in low frequency 
radio observations. In three earlier papers ( Bharadwaj, Nath & Sethi, 2001; 
Bharadwaj & Sethi, 2002; and Bharadwaj & Pandey, 2003; hereafter referred to 
as Paper a, b and c respectively), and in the present paper we investigate issues 
related to calculating the expected signal and detecting it. 

We propose (Papers b and c ) that the optimal observational strategy for 
detecting this signal is to deal directly with the complex visibilities measured 
in radio interferometric observations. Briefiy introducing the terms involved, 
we remind the reader that the quantity measured in radio interferometric ob- 
servations with an array of antennas is the complex visibility V{\], v). This is 
measured for every pair of antennas at every frequency channel in the observa- 
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tion band. For every pair of antennas it is convenient to express the visibility as 
a function of U = d/A ie. the separation between the two antennas d expressed 
in units of the wavelength A. We refer to the different possible values of U as 
baselines. One of the big advantages of dealing directly with the visibilities is 
that the system noise contribution to the visibilities is uncorrelated. The visi- 
bilities respond only to the fluctuations in the redshifted HI emission. In Paper 
b we showed that the correlation expected between the visibilities ^(U, v) and 
V{IJ, v+Au) measured at the same baseline at two slightly different frequencies 
is 



where k = ^(2 tt J7/r)2 + fc^, Iq is the specific intensity expected from the 
HI emission if the HI were uniformly distributed, 6q = 0.6 x ^fwhm, ^^fwhm 
being the angular width of the primary beam of the individual antennas, r 
is the comoving distance to the HI from which the radiation originated, and 



tribution in rcdshift space at the epoch when the HI emission originated. 

To summarize, the visibility-visibility cross-correlation (hereafter rcfcrccd to 
as the visibility correlation) directly probes the power spectrum of HI fluctua- 
tions at the epoch where the HI emission originated. This holds the possibility 
of allowing us to study the large scale structures at high rcdshifts. A point 
to note is that the the visibility correlations at a baseline U receives contribu- 
tion from the power spectrum only for Fourier mode k > fcinin = (27r/r)C/, and 
for the CDM-likc power spectrum most of the contribution comes from Fourier 
modes around fcmin. So, it may be said that the correlations at a baseline U 
probes the power spectrum at the Fourier mode (27r/r)C/. 

In the earlier work we treated the HI as being continuously distributed 
whereas in reality the HI resides in discrete gas clouds. In addition, it was as- 
sumed that the HI distribution is an unbiased representation of the underlying 
dark matter distribution, and we used the linear theory of density perturbations 





PD^P{k) l + 



is the power spectrum of the fluctuations in the HI dis- 
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to follow the evolution of fluctuations in the dark matter distribution. It is these 
assumptions which allow us to express the power spectrum of fluctuations in the 
HI distribution in redshift space at the epoch when the HI emission originated 
(in eq. 1) in terms of b the linear bias parameter (taken to be 1), £) the grow- 
ing mode of linear density perturbations (Peebles, 1980) at the epoch when 

the HI emission originated, P{k) the present power spectrum of dark matter 

r 

k 

density fluctuations calculated using linear theory and the factor 1 + P-^ 
which takes the power spectrum from real space to redshift space in the linear 
theory of redshift distortions (Kaiser 1987) . Here we report progress on two 
counts. First, we have used a PM N-body code to evolve the fluctuations in 
the dark matter distribution, thereby incorporating possible non-linear effects. 
Second, we have assigned HI masses to the dark matter particles in the N-body 
code and this was used to simulate the redshifted HI emission. So we have also 
been able to incorporate the fact that the HI gas is contained in discrete clouds. 
The predictions for the HI signal expected at GMRT presented in this paper 
incorporate both these effects. We still retain the assumption that the HI is an 
unbiased tracer of the dark matter. 

We next present a brief outline of this paper. In section 2 we discuss the 
method that was used to simulate the HI signal, and in Section 3. we present 
the results of our investigations. In Section 4 we discuss the results and present 
conclusions. 

Finally, it should be point out that there been alternative lines of approach 
investigating the possibility of using HI observations to study large scale struc- 
tures at 2 ~ 3 (Sunyaev & Zeldovich, 1975; Subramanian & Padmanabhan, 
1993; Kumar, Padmanabhan & Subramanian, 1995; Weinberg et al., 1996; 
Bagla, Nath & Padmanabhan, 1997; Bagla & White, 2002). The reader is 
referred to Papers a and b for a detailed comparison of these approaches with 
that adopted here. 
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Table 1. 



z/c (MHz) 


z 


^FWHM 


r (Mpc) 


L (Mpc) 


AL (Mpc) 


-/Vdm 






610 


1.33 


0.9° 


4030 


128 


0.5 


128^ 


4 


19 


325 


3.37 


1.8° 


6686 


512 


1 


256^ 


4 


9 



2 Methodology. 

We have simulated the visibility correlations expected at two of the GMRT 
observing frequency bands centered at = 610 MHz and 325 MHz. The simu- 
lations were carried out in three steps 

1. Using a PM N-body code to simulate the dark matter distribution at the 
redshift where the HI emission originated 

2. Assigning HI masses to the particles used in the N-body code and calcu- 
lating the flux expected from each HI cloud 

3. Calculating the complex visibilities arising from the distribution of HI 
clouds and computing the visibility correlations. 

We next discuss the salient features of each of these steps. The values h = 
0.7, fimo = 0.3 and Qao = 0.7 were used throughout. 

2.1 The N-body Simulations. 

We have used a Particle-Mesh (PM) N-body code to simulate the dark matter 
distribution at the redshift z where the HI emission originated. The simulation 
volume was a cubic box of comoving volume li^ . The size L was chosen so that 
it is approximately twice the comoving distance subtended by ^fwhm of the 
GMRT primary beam. 

The values of r the comoving distance to the region from where the HI 
emission originated, the grid spacing of the mesh AL, and the number of dark 
matter particles used in each simulation A^dm are all shown in Table 1. 
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The initial power spectrum of dark matter density fluctuations at Zin (shown 
in Table 1) is normalized to COBE (Bunn & White 1996), and its shape is 
determined using the analytic fitting form for the CDM power spectrum given 
by Efstathiou, Bond and White (1992). The value of the shape parameter turns 
out to be r = 0.2 for the set of cosmological parameters used here. We have 
run the N-body code for A^sim (Table 1) independent realisations of the initial 
conditions and the final results for the visibility correlations were averaged over 
all the realisations. 
10000 
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Figure 1: This shows the power spectrum of density fluctuations in the dark 
matter distribution a,t z = 1.33. The vertical lines show the smallest Fourier 
mode fcmin = {2iT/r)U which contributes to the visibility correlations at a base- 
line U. This is shown for the different values of U indicated in the figure. 




The N-body code gives the final positions and peculiar velocities of the A^dm 
dark matter particles in the simulation. The power spectrum of the density fluc- 
tuations in the dark matter distribution a,t z = 1.33 and z = 3.37 are shown in 
Figures 1 and 2 respectively. In both the flgures we have shown the power spee- 
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Figure 2: This is the same as figure 2 except that it is at z = 3.37 which 
corresponds to 325MHz. 

trum for the range of Fourier modes which will make a significant contribution 
to the visibility correlation at the baselines where the signal is expected to be 
strongest. We find that a,t z = 1.33 (610MHz) the power spectrum obtained from 
the N-body simulation shows substantial differences from the power spectrum 
calculated using linear theory at Fourier modes k > 0.3Mpc~^. Converting 
to baselines, we expect non-linear effects to be important for U > 200. At 
z = 3.37 (325MHz) there are differences between the N-body and linear power 
spectrum at fc > 0.3Mpc~^, but the non-linear effects are not as pronounced as 
at 325MHz. Converting to baselines, non-linear effects will influence the signal 
at baselines U > 400. 

2.2 Assigning HI masses. 

We model the HI clouds as randomly oriented rotating disks of radius R, column 
density Nui and rotation velocity V. The values of R and V are held flxed 
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Table 2. 



(MHz) 


a 


R Kpc 


n^i(z)Mpc-^ 


^clouds 


610 


1.2 


10 


2.0 X 10"^ 


42467 


610 


1.2 


8 


3.2 X 10-2 


66354 


610 


1.2 


5 


8.1 X 10-2 


169865 


610 


1.2 


2 


5.1 X 10-1 


1061651 


610 


1.7 


10 


3.7 X 10-2 


76764 


610 


1.7 


8 


5.7 X 10-2 


119943 


610 


1.7 


5 


1.5 X 10-1 


307055 


610 


1.7 


2 


9.2 X 10-1 


1919088 


325 


0.8 


10 


1.3 X 10-2 


1799234 


325 


1.2 


10 


2.0 X 10-2 


2717825 


325 


1.2 


8 


3.2 X 10-2 


4246602 


325 


1.2 


5 


8.1 X 10-2 


10871300 



in each simulation, and we have run simulations with different sets of values 
for these parameters. It is assumed that the column densities have a power 
law distribution in the range 2 x 102" < iVHi/(atoms/cm2) < 1 x IO22 and the 
comoving number density of HI clouds with column densities in the interval rfA^ni 
is BN^" dNiii- The total comoving number density of HI clouds at redshift z 
is 

/■AfHl[max] 

n'k,{z)=B / N^rdNm (2) 

J Nm [min] 

and the comoving mass density is 

/•AfHl[max] 

p'hi{z) = B / {TTR^Nmmm)N^rdNm (3) 

<y N-ni [min] 

where TOhi is the mass of the hydrogen atom. The normalisation coefficient B 
is determined by using equation (3) to calculate i}gas{z) 



^gas{z) = ^flm{z) = -^Pm{z) (4) 
We use ilnas = in all the simulations. 



This model fixes the total number of HI clouds in the simulation volume 
■^clouds = L^n-iiiiz) . The total number of clouds scales as 

1 — a 

A^clouds OC i?-2 (5) 

2 — a 
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as we vary a the slope of the column density distribution or the radius of the 
clouds R. For large values of R the HI is distributed in a few clouds with large 
masses, whereas there are many clouds with low HI masses when R is small. Our 
model has three free parameters, namely a, R and V. We have run simulations 
varying a and R (Table 2) for V = 100 km/sand 200 km/s. 

We randomly select A^ciouds particles from the output of the N-body sim- 
ulation and these are identified as HI clouds. The HI mass of each cloud is 
-^Hi = T^R^Nuimm, where the column density is drawn randomly from the 
power-law distribution discussed earlier. The center of the simulation volume 
is aligned with the center of the GMRT primary beam and it is located at a 
comoving distance corresponding to the redshift z. The comoving distance to 
each cloud is used to calculate its angular position and redshift. The redshift 
is used to determine the luminosity distance which is used to calculate the flux 
from the individual clouds. The effect of the peculiar velocity is incorporated 
when calculating Vg the frequency at which the HI emission from each cloud is 
received. The line width Av of the HI emission line from each cloud is calcu- 
lated using Av =1 sin^ | 2i'oV/c, where 2uoV/c is the line width if the disk of 
the galaxy were viewed edge on and 6 is the angle between the normal to the 
disk and the line of sight. 

To summarize, at the end of this stage of the simulation we have -A^ciouds HI 
clouds. For each cloud we have its angular position 0°',and the flux density F", 
frequency and line-width Ai/" of the redshifted HI emission. Here the index 
a {1 < a < A^ciouds) refers to the different clouds in the simulation. 

2.3 Calculating visibility correlations 

We flrst describe how we have calculated the complex visibilities that would 
be measured in GMRT radio observations of the HI distribution generated in 
the simulation. The observations are carried out at NC frequency channels 
{z/i, 1/2, f3, z^Afc} covering a frequency band B centered at the frequency Vc- 
We have used B = 8MHz and NC = 64 at i^c = 610MHz, and B = 8MHz and 
NC = 128 at i/c = 325MHz. 
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For the purpose of this paper we assume that the antennas are distributed 
on a plane, and that they all point vertically up wards. The beam pattern 
A{9) quantifies how the individual antenna, pointing up wards, responds to 
signals from different directions in the sky. This is assumed to be a Gaussian 
A{e) = e-^'/^S where 60 = 0.6 x 6Ifwhm (Table 1). 

The position of each antenna can be denoted by a two dimensional vector dj. 
The quantity measured in interferometric observations is the visibility ^(U, v) 
which is recorded for every independent pair of antennas (baseline) at every 
frequency channel in the band. For any pair of antennas, the visibility depends 
on the vector d = dj — dj joining the position of the two antennas. It is 
convenient to express the visibility as a function of the variable U which is d 
expressed in units of the wavelength i.e. U = d/A. The signal arising from the 
clustering pattern of the HI clouds will be strongest at the small baselines, and 
our calculations have been limited to this. We have considered a square grid 
of baselines extending from —Umax to f/max with resolution 6U. We have used 
Umax = 400 and SU = 10. The complex visibility has been calculated for each 
baseline U on the grid using 

V{U,p)= ^ A{r ) F«e-'^'^^-^"'°0( ' ''^ J° ' ) (6) 

a=l 

where the function 0{x) is defined such that 0{x) = 1 for a; < 1, else 0{x) = 0. 
It is to be noted that in an actual GMRT observation the baselines will have a 
complicated distribution depending on which part of the sky is observed and the 
duration of the observation. Given the fact that the signal we axe interested in 
is statistical in nature, and that we are interested in making generic predictions 
about the signal expected in a typical GMRT observation, a square grid of 
baselines is adequate. 

The final step in the simulation is to calculate the visibility correlation 
{V{U,i')V*(U,i/ + Av)). The angular brackets () indicate the ensemble av- 
erage, and we have averaged over the A^sim different realisation of the N-body 
simulation. In addition, the correlation depends only on the separation in fre- 
quency I Av I, and the magnitude U =\ IJ \. So, for a fixed values of Av and 
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U we have averaged over all possible pairs of frequencies and baselines which 
match these values. 

The analytic calculations (Papers b and c) where the HI is assumed to have 
a continuous distribution, predict the imaginary part of the visibility correlation 
function to be zero, and the clustering signal is manifest in only the real part. In 
the simulations we get a very small, but non-zero imaginary component. This 
is not discussed in the rest of the paper where we present results for the real 
component only. 

3 Results 

In this section we present results for the visibility correlation as obtained from 
our simulations. We compare these with the analj^tic predictions of Papers b 
and c and investigate the effect of two factors (1.) the non-linear evolution of 
the density fluctuations, and (2.) the discrete nature of the HI distribution. To 
get a better understanding of the second effect, we present results varying the 
parameters of the HI distribution. 

3.1 610 MHz 

Figure 3 shows the visibility correlations for U = 100, the results at smaller 
baselines show a similar behaviour. The visibility correlation at the baseline U = 
100 receives contributions mainly from Fourier modes around k ~ 0.2Mpc~^ 
which is in the linear regime (Figure 1), and we expect a good agreement with 
the analytic, linear predictions. We find that for Au < 0.5MHz the analytic 
predictions are larger than the correlations obtained in the simulations, and this 
is reversed for Av > 0.5MHz. This discrepancy can be explained if we take into 
account the fact that visibility correlation actually responds to the clustering 
in redshift space. It is known (e.g. Suto & Suginohara 1991; Graman, Gen & 
Bahcall, 1993; Fisher et al. 1994; Brainerd et al., 1996; Bromley, Warren & 
Zurek, 1997) that non-linear effects can be important in redshift space even on 
scales where the clustering in real space is well described by linear perturbation 
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Figure 3: This shows the correlation expected between the visibihties ¥{13,1/) 
and 1^(U, u + Av) at the same basehnes at two different frequencies. The 
rotational velocity of the HI disk is assumed to be y = 200km/s. The other 
parameters of the HI distribution take on values shown in the figure. These 
results are for the 610MHz band. 

theory. It has been shown that this can be modeled by taking into account 
the effect of the random motions along the line of sight (e.g. Fisher et aL1994, 
Peacock & Dodds 1994, Ballinger, Peacock & Heavens 1996). We incorporate 
this by multiplying the power spectrum with exp[— fc^cr^] in our analytic formulas 
for the visibility correlation. This gives the modified formula 



for the visibility correlation. We find that for a = 200km/s/i?o this gives a 
good fit to the results of the simulations, and this is also shown in the figure. 
We next shift our attention to how the results depend on the parameters of 
the HI distribution. We find that for r = lOKpc where the bulk of the HI 





(7) 
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is distributed in a few clouds with large HI masses the results show a 20% 
increment at small values of l\v compared to the models with smaller values of 
r. This excess correlation at small Az/ arises from the fact that the HI emission 
from an individual cloud will be spread across a width 8v in frequency. The 
correlation between the HI emission from the same HI cloud at two different 
frequencies will contribute to the visibility correlations when Ai^ < bv. The 
contribution from this signal is significant in comparison to that arising from 
the clustering of the HI clouds when the total HI is distributed in a few clouds 
with large HI masses each. The contribution to the visibility correlation from 
within individual HI clouds goes down as r is reduced and the HI is distributed 
among many clouds each with small HI masses. There is very little difference 
between the results for r = 5Kpc and 2Kpc and we may treat this as the result 
if the HI were continuously distributed. 

The results in Figure 3 are for the rotational velocity V = 200 km/s. We have 
also done simulations using V — 100 km/s. We find that there are differences 
(< 20%) only at small values of Ai/. The effect of decreasing V is to decrease 



13 



the frequency width of the HI emission hne from individual clouds which results 
in a higher value of the HI flux density. This does not effect the clustering 
signal but enhances the contribution to the visibility correlation arising from 
the emission of a single HI cloud. As changing V does not affect the results very 
much, in this subsection we show the results for V = 200km/s only. 

Figure 4 shows the results for U = 300. We find that the discrepancy be- 
tween the linear, analytic predictions and the results of our simulations increases 
at larger values of U. Except at very small values of Av, the simulated values 
are larger than the linear predictions. This is because the larger baselines probe 
smaller length scales which are significantly nonlinear (figure 1) and the ampli- 
tude of the fluctuations is larger than predicted by linear theory. At smaller 
scales the fluctuations are non-linear even in real space, and the modified formula 
(equation 7) based on only redshift space considerations grossly underestimates 
the visibility correlations. An important point is that at larger values of U the 
visibility correlations calculated in the simulations do not fall as sharply with 
increasing Av as predicted in the linear calculations. In Paper c we found that 
the visibility correlations decay as oc expl—Av/K], where the decay constant 
varies as if oc U~^'^ i.e. the decay is faster at larger baselines. Our simula- 
tions show that the decay with increasing Au is slower than predicted using 
linear theory. This is a consequence of the fact that the density fluctuations are 
non-linear on the length-scales being probed at these baselines. Another point 
to note is that the dependence on the parameters R and a, becomes relatively 
more pronounced at large values of U. 

3.2 325 MHz 

A point which should be mentioned right at the start is that at 325MHz we 
are restricted in the values of R for which we are able to carry out simulations. 
At 325MHz the simulation volume is pretty large (Table 1) and for R = 2Kpc 
the total number of HI clouds in the simulation volume becomes too large for 
our computational resources. Also, in this subsection we use V = lOOkm/s 
in our simulations. Figure 5 shows the results for the visibility correlations at 
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2.5e-09 



325 MHz 
U=100 




-5e-10 



Linear 
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1.5 



AV Mhz 



Figure 5: This shows the expected correlation between the visibilities ¥{13,1/) 
and 1^(U, u + Av) at the same baselines at two different frequencies. The 
rotational velocity of the HI disk is assumed to be y = lOOkm/s. The other 
parameters of the HI distribution take on values shown in the figure. These 
results are for the 325MHz band. 



U = 100. The behaviour at smaller baseline'es is not very different. We find that 
for Ai/ < IMHz the predictions of the analytic, linear calculations (Paper c) are 
very close to the values obtained in the simulation for r = 8Kpc and a = 1.2. 
The results of the simulation are slightly larger than the analytic predictions 
when r = lOKpc, and they are somewhat smaller than the analytic predictions 
at r = 5Kpc. For Av > IMHz the results of the simulation are the same for all 
the parameters, and the value is slightly more than the analytic prediction. The 
power spectrum (Figure 2) is in the linear regime at the Fourier modes which 
contribute to the visibility correlations at U = 100. The discrepancy between 
the analj^ic predictions and the results of our simulations can be attributed 
to a combination of the two factors discussed earlier (1.) the effect of random 
motions on the redshift space clustering and (2.) correlations between the HI 
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Figure 6: This is the same as Figure 5 except that it shows results for U = 400. 



emission from the same cloud at different frequency channels. 

Figure 6 shows the visibility correlations at U = 400. We find that the 
behaviour of the visibility correlations docs not change very much for baselines 
in the range 100 < U < 400. The power spectrum (Figure 2) starts getting 
non-linear at Fourier modes corresponding to U = 400, but the effect is not 
very significant. 

4 Discussion and Conclusions. 

We take up for discussion two issues pertaining to the way we have modeled the 
distribution of HI clouds. First is our assumption that the HI clouds responsible 
for damped Lyman-a absorption lines are rotating disks, all with the same 
radius and rotational velocity. Prochaska and Wolfe (1998) have proposed that 
these HI clouds are the gaseous progenitors of present day galaxies. They have 
attempted to explain the observed kinematic of the damped Lyman-a absorption 
lines using a thick, rotating disk model for the HI clouds. Another model ( 



16 



Haehnelt, Steinmetz &Rauch, 1998) proposes that the observations could be 
better explained by modeling the absorption systems as protogalactic gas clumps 
undergoing merger. The second issue is our assumption that the HI column 
densities have a power law distribution and our choice of the value of the index 
a. Lanzetta et al. (1991) show that the column density distribution at 2; ~ 2.5 
can be described by a power law with a ~ 1.7. In a later paper Lanzetta, Wolfe 
and Turnshek (1995) show that the column density distribution evolves quite 
strongly with redshift, there being a tendency toward more high column density 
clouds at higher redshifts. In our work we have run simulations for two values 
of the index i.e. a = 1.2 and 1.7. Having very briefly reviewed some of the 
prevalent views and having compared our assumptions with them, we note that 
our simulations seem to indicate that the visibility correlation signal does not 
depend very critically on the details of the properties of the HI clouds. We find 
that the visibiHty correlation signal has contribution from mainly two effects (1.) 
correlations caused by the emission from the same cloud to different frequency 
channels, and (2.) the clustering of the clouds in redshift space. The first effect 
is seen in the correlation at small values of Au, where Av is smaller than the 
width of the HI line from an individual cloud. This effect manifests itself as a 
rise in the visibility correlations at small values of Az/(< 0.5MHz). This effect 
is enhanced if the HI is distributed in a few clouds with large HI masses as 
compared to the situation where the HI is in many small clouds with low HI 
masses. This seems to be the only effect of the fact that the HI is distributed 
in discrete HI clouds and is not continuously distributed. This is also the only 
place where the details of the HI distribution affects the visibiHty correlation. 

We next turn our attention to the contribution to the visibility correlation 
signal from the clustering of the HI clouds. We may take the results for the 
visibility correlations at the values of the parameters R and a where A^ciouds 
is maximum as representing the results when the discrete nature of the HI 
distribution can be neglected (continuum limit). This assumption is justified at 
610 MHz where the is very little difference in the results between R = 5 Kpc and 
R = 2 Kpc. We have not carried out simulations for ii < 5 Kpc at 325 MHz, and 
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simulations with smaller values of R are needed at this frequency before we can 
be sure that the results for i? = 5 Kpc really represent the continuum limit. Let 
us first discuss the results at the baselines for which the visibility correlation 
probes the power spectrum at length-scales which are in the linear regime. This 
is true for baselines with U < 100 at 610 MHz (Figure 1). At 325 MHz most of 
the baselines which we have studied probes the power spectrum in the linear 
regime (Figure 2). For all these baselines we find that the simulated values are 
less than the predictions of linear theory at small Av and the simulated values 
are larger than the linear predictions at large Au. The transition occurs in the 
range Av ~ 0.5 — 1 MHz. We propose that this discrepancy is a consequence 
of the fact that the fluctuations in the HI distribution in redshift space may be 
non-linear even on length-scales where linear theory holds in real space. This 
can be modeled by incorporating the effect of random peculiar velocities on 
the redshift space HI distribution. We show that including this effect gives a 
good fit to the simulated results atU = 100 for 610 MHz. At larger baselines the 
visibility correlation probes the power spectrum on length-scales where it is non- 
linear. Non-linear effects start influencing the visibility correlation at baselines 
U > 200 for 610 MHz, and these effects are very significant hy U = 400. As a 
consequence of these effects the simulated visibility correlations do not fall of 
with increasing Av as quickly as predicted by linear theory. Also, the simulated 
values are larger than the linear predictions everywhere except at very small 
values of Av. The range of Av where the simulated value are less than the 
linear predictions decreases with increasing U. 

In conclusion we note that the HI signal predicted by our simulations are 
not drastically different from the analytic predictions presented earlier. In this 
paper we have been able to address the effects of the discrete nature of the HI 
distribution and the non-linear nature of the HI fluctuations in redshift space. 
We now have the tools necessary to simulate the HI signal expected at the 
GMRT. A full simulation of a GMRT observation requires us to also include 
the system noise as well as various galactic and extragalactic radio sources. 
Only then will we be able to make definite predictions as to whether it will be 
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possible to detect the HI signal or not. Work is correctly underway on this. The 
preliminary results indicate that it will be posiible to have a. 5 a detection at 
610 MHz with one thousand hours of observation. 
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